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ABSTRACT 

Star formation theory predicts that short-period M-dwarf binaries with highly 
unequal-mass components are rare. Firstly, the mass ratio of close binary systems 
is driven to unity due to the secondary preferentially accreting gas with high angu- 
lar momentum. Secondly, both dynamical decay of multiple systems and interactions 
with tertiary stars that tighten the binary orbit will eject the lowest mass member. 
Generally, only the two most massive stars are paired after such interactions, and the 
frequency of tight unequal-mass binaries is expected to decrease steeply with primary 
mass. In this paper we present the discovery of a highly unequal-mass eclipsing M- 
dwarf binary, providing a unique constraint on binary star formation theory and on 
evolutionary models for low-mass binary stars. The binary is discovered using high- 
precision infrared light curves from the WFCAM Transit Survey (WTS) and has an 
orbital period of 2.44 d. We find stellar masses of Mi = 0.53 ± 0.02 M and M 2 = 
0.143 ± 0.006 M Q (mass ratio 0.27), and radii of R x = 0.51 ± 0.01 R Q and R 2 = 0.174 
± 0.006 Rq. This puts the companion in a very sparsely sampled and important late 
M-dwarf mass-regime. Since both stars share the same age and metallicity and strad- 
dle the theoretical boundary between fully and partially convective stellar interiors, a 
comparison can be made to model predictions over a large range of M-dwarf masses 
using the same model isochrone. Both stars appear to have a slightly inflated radius 
compared to 1 Gyr model predictions for their masses, but future work is needed to 
properly account for the effects of star spots on the light curve solution. A significant, 
subsynchronous, ^2.56 d signal with ~2% peak-to-peak amplitude is detected in the 
WFCAM light curve, which we attribute to rotational modulation of cool star spots. 
We propose that the subsynchronous rotation is either due to a stable star-spot com- 
plex at high latitude on the (magnetically active) primary (i.e. differential rotation), 
or to additional magnetic braking, or to interaction of the binary with a third body 
or circumbinary disk during its pre-main-sequence phase. 

Key words: stars: low mass - stars: starspots - binaries: eclipsing 



1 INTRODUCTION 

Stellar evolution models find it difficult to accurately 
reproduce the fundamental properties of M-dwarf stars 
(M* <O.65M0), which are the most abundant population 
of stars in the Milky Way (over 70% in number; Henry 
et al. 1997). Furthermore, the formation process of binary 
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M-dwarfs and their migration to close orbits are not well- 
understood (e.g. Goodwin & Whitworth 2007; Clarke 2012; 
Nefs et al. 2012). M-dwarfs are important to astrophysics be- 
cause they help us understand a variety of problems, from 
the local star formation history, to the shape of the stellar 
(initial) mass function (e.g. Reid 1999; Gizis, Reid & Hawley 
2002; Bochanski et al. 2007). Moreover, M-dwarfs are now 
being recognised as prime targets in the hunt for Earth-like 
exoplanets (e.g. Martin et al. 2005; Nutzman & Charbon- 
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neau 2008), as their small stellar size results in deeper tran- 
sits and their smaller stellar mass causes larger reflex mo- 
tions induced by planetary companions. However, the cur- 
rent uncertainty in our understanding of M-dwarf formation 
and evolution means that the parameters of their planets, 
which scale with those of the host star, can not be deter- 
mined to high accuracy, limiting a detailed characterisation 
of their compositions and atmospheres. 

Hydrodynamical simulations (e.g. Bate, Bonnell & 
Bromm 2002) show that when a star-forming cloud col- 
lapses, it fragments into low-number multiple systems, which 
can be broken up by dynamical interactions (e.g. Goodwin 
et al. 2007). It is predicted that short-period M-dwarf bina- 
ries with significantly unequal-mass components should be 
very rare. Firstly, this hypothesis is motivated by the Bate, 
Bonnell & Bromm (2002) simulations that suggest that in- 
falling gas with high angular momentum is preferentially 
accreted onto the lower-mass component of a binary during 
formation, driving the mass-ratio to unity (see also Bate & 
Bonnell 1997; Bate 2000). Close binaries in particular are 
expected to experience significant accretion of gas in or- 
der to migrate to tighter orbits. The angular momentum 
of the gas increases as the accretion proceeds and is higher 
for closer binary systems, which accrete more gas relative to 
their initial masses than wider binaries. Whereas gas with 
low angular momentum is mainly accreted by the primary 
because it essentially falls straight onto the centre of mass of 
an unequal-mass system, high angular momentum gas falls 
in further away, closer to the orbital radius of the secondary, 
and need not to gain as much momentum to be captured by 
the companion (e.g. Bate & Bonnell 1997). 

An enhanced preference of equal mass ratios is expected 
especially for tight low-mass binaries. This is because dy- 
namical decay of multiple systems and exchange interactions 
with single stars in the collapsing cloud, which increase the 
binary mass ratio and at the same time tighten the binary 
orbit, are biased against the lowest mass stars. In dynam- 
ical decay typically the least massive component is ejected 
on a short timescale, due to the instability of multiple star 
systems (e.g. Anosova 1986). Furthermore, in dynamic ex- 
change interactions the lowest mass star is removed and re- 
placed by the higher mass intruder. This means that only 
the two most massive stars survive an interaction, indicating 
that both the occurrence frequency of binary and multiple 
systems and the number of close and unequal-mass binaries 
steeply decrease with decreasing primary mass. Therefore 
each of the mechanisms involved in producing close and low- 
mass binaries favors the production of equal-mass systems 
and highly unequal, short-period, M-dwarf binary systems 
should thus be rare. 

There is observational evidence for this hypothesis. Wis- 
niewski et al. (2012) recently propose a lack of unequal-mass 
stellar binaries at periods shorter than ~100 d, combin- 
ing current results from radial velocity, transit and imag- 
ing studies. Delfosse et al. (2004) find that M-dwarf binaries 
with orbital period P or b < 50 d possess a mass-ratio dis- 
tribution which is peaked around 1 ('twins'), whereas wider 
binaries have a flat distribution. Clarke, Blake & Knapp 
(2012) find a frequency of close M-dwarf binary stars with a 
separation of less than 0.4 AU of only 3.0^q g%, and argue 
that the frequency of <0.4 AU binaries is decreasing from 
26% at 1OM to 1% at 0.1M o (see e.g. Lada (2006) and 



Raghavan et al. 2010). Their findings suggest that low-mass 
stars not only have lower multiplicity but their tightest sys- 
tems are also intrinsically rare. Bouchy et al. (2011) propose 
that G-type or lower mass stars have stronger disk braking 
than more massive stars, which causes any short-period low- 
mass companion to migrate inwards and become engulfed by 
the primary. 

Detached, double-lined eclipsing binaries (EBs) provide 
a model-independent method for accurately calibrating the 
formation and evolution of stars (Andersen et al. 1991; Tor- 
res, Andersen & Gimenez 2010). Dynamical measurements 
of the masses and radii and empirical measurements of the 
temperatures of M-dwarf stars in eclipsing binary systems, 
suggest that current models under-predict the radii of M- 
dwarfs by 5 — 15%, and over-predict their effective temper- 
atures by 3 — 5% (typically 100-200K; e.g. Lopez-Morales 
& Ribas 2005; Lopez-Morales 2007; Torres, Andersen & 
Gimenez 2010). These discrepancies are attributed to pos- 
sible metallicity variations (e.g. Berger et al. 2006), but 
more likely they are due to strong magnetic fields and spots 
present on the stars of the observed tidally-locked, short- 
period M-dwarf binaries (e.g. Mullan & MacDonald 2001; 
Ribas 2006; Chabrier, Gallardo & Baraffe 2007). However, 
there is no current model that can accurately reproduce all 
of the observed dynamical measurements. 

It is thought that stars with masses M < 0.35Mq 
converge towards agreement with the current models be- 
cause they are likely to have fully-convective atmospheres, 
and thus suffer less from the inflating effects of magnetic 
inhibition (Kraus et al. 2011). However, even the longest- 
period, non-synchronised M-dwarf EBs (MEBs) with sec- 
ondary components in the fully-convective regime (Irwin et 
al. 2011; Doyle et al. 2011), still show significant radius 
inflation, despite much lower magnetic fields (see Birkby 
et al. 2012). For the lowest mass main-sequence M-dwarfs 
(0.08Mq < M < O.2M0), there is even more uncertainty, 
due to a paucity of model-independent dynamical measure- 
ments. Only a few young objects in the Orion Nebula have 
data (see Irwin et al. 2007; Stassun, Mathieu & Valenti 
2007) . Constraints on M-dwarf evolution isochrones are fur- 
ther hindered by the apparent preference for equal mass 
MEBs. Low-mass ratio MEBs are valuable because their 
shared age and metallicity allow a more stringent assess- 
ment of the stellar model predictions over the wide span of 
M-dwarf masses. 

In this paper we describe the discovery and charac- 
terisation of a main-sequence unequal mass, short-period, 
detached M-dwarf eclipsing binary system, whose compo- 
nents straddle the fully-convective boundary, and whose sec- 
ondary star resides in the uncharted < 0.2Mq mass regime. 
The binary is discovered using high-precision infrared light 
curves from the WFCAM Transit Survey (WTS; Birkby, 
Hodgkin & Pinfield 2011; Birkby et al. 2012; Kovacs et al., 
submitted). The WTS is an ongoing photometric monitor- 
ing campaign that operates as a back-up program running 
on the 3.8m United Kingdom Infrared Telescope (UKIRT) 
at Mauna Kea, Hawaii. By observing in the infrared, the 
WTS is optimised for precision photometry of cool low-mass 
stars. Its primary objective is to hunt for planets orbiting 
M-dwarfs by regular monitoring of ~6000 early- to mid M- 
dwarfs (for J<16), but also to characterise the host stars. 

In Section 2 we present the observations and the data 
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reduction of infrared and optical time-series of the eclipses of 
our binary, WTS 19g-4-02069, and present low- to medium 
resolution spectroscopy in the optical and in the infrared ri- 
band. In Sections 3 and 4 we characterise the components 
of WTS 19g-4-02069 using the available data, obtaining in- 
dividual masses, radii, effective temperatures and activity 
indicators. In Section 5 we discuss the significance of the 
binary in the context of current low-mass stellar evolution 
models. 



2 OBSERVATIONS & DATA REDUCTION 
2.1 WTS J-band time series photometry 

The WTS, in operation since 2007 August, is awarded 200 
nights of observing time on UKIRT. The survey uses the 
UKIRT Wide Field Camera (WFCAM), which has four 
2048x2048 18 fim HgCdTe Rockwell Hawaii-II imaging ar- 
rays that each cover 13.65x13.65 arcmin 2 on sky (with a 
pixel resolution of 0.4 arcsec pixel -1 ), and are seperated by 
94 per cent of a chip width (Casali et al. 2007). Observa- 
tions for the WTS are obtained in the J-band (1.25 /xm), 
which maximises the sensitivity to M-dwarfs with effective 
temperatures T e ff <4000K. 

The survey targets four 1.5 deg 2 fields centered around 
3h, 7h, 17h and 19h in Right Ascension, selected to give 
both year-round visibility, an optimal number of dwarfs ver- 
sus giants, relatively low reddening [E(B-V) between 0.057 
and 0.234] and reduced contamination by blending stars, 
by observing close to but outside the galactic plane (galac- 
tic latitude b >5°). The infrared light curves have an aver- 
age cadence of 15 min. For each field, single deep exposures 
in the full WFCAM ZYJHK system are also obtained, to 
aid the photometric identification of M-dwarfs through fit- 
ting of the broad band spectral energy distributions (SEDs). 
The SED for WTS 19g-4-02069, extending from Sloan Dig- 
itized Sky Survey (SDSS) u band to infrared Wide Field 
Infrared Survey Explorer (WISE) 4.6/im, is shown in Table 
1. The observing strategy, pipeline data reduction and WTS 
light curve generation is extensively described in Birkby et 
al. (2012) and Kovacs et al., submitted and the interested 
reader is referred to these two publications for more details. 

The subject of this paper is selected from the list of 
16 well-sampled detached M-dwarf eclipsing binaries with 
J sj 16 as presented in Birkby et al. (2012) in the 19hr field, 
the WTS target field which currently has the most extensive 
observational coverage (~1200 epochs). An initial source 
detection is performed using the Box-Least-Squares (BLS) 
algorithm, OCCFIT (Aigrain & Irwin 2004; Miller et al. 
2008). Fitting of the broad band SED of WTS 19g-4-02069 
yields a system effective temperature of T e // ~3050K, the 
lowest of the Birkby et al. (2012) sample, indicating a low- 
mass eclipsing M-dwarf binary system. The out-of-eclipse 
root mean square (rms) scatter of WTS 19g-4-02069 is rel- 
atively high compared to other stars of similar magnitude 
(~12.5mmag per datapoint), whereas ~6mmag is expected. 
We attribute this to stellar activity (see the discussion in 
Section 4.1.1). The WTS J band data for WTS 19g-4-02069 
are given in Table 2 and shown in Figure [1] folded on the 
binary orbital period. 



OJ2000 19:35:03.55 
5j2000 +36:31:16.49 
Broad-band SED WTS 19g-4-02069 



Band name 


Central A(/im) 


magnitude 


u 1 


0.35 


23.594(±3.0) 


g 


0.47 


20.262(±0.02) 


r 


0.62 


18.866(±0.01) 


i 


0.75 


17.277(±0.01) 


z 


0.89 


16.335(±0.01) 


Z 


0.89 


15.928(±0.005) 


Y 


1.03 


15.426(±0.007) 


J 


1.25 


14.843(±0.004) 


11 


1.63 


14.271(±0.004) 


K 


2.20 


13.952(±0.004) 


WISE1 


3.4 


13.842(±0.028) 


WISE2 


4.6 


13.812(±0.044) 



Table 1. Broad-band spectral energy distribution for the eclips- 
ing M-dwarf binary WTS 19g-4-02069. SDSS u,g,r,i and z magni- 
tudes are quoted in the AB magnitude system, whereas the WF- 
CAM Z,Y,J.H and K magnitudes are in the Vega system. X SDSS 
u is uncertain because of a red leak (Abazajian et al. 2004). The 
entries WISE1 and WISE2 refer to the first two wavelength chan- 
nels of the WISE at 3.4 and 4.6/xm. The source is too faint for 
detection at 12 and 22/xm. 



WFCAM HJD J WTS ctJwts 

J-band -2454000 

(days) (mag) (mag) 

317.808593 14.8003 0.0050 

INT HJD AmnpfT crmnpfT 

i-band (mag) (mag) 
Primary event 

1403.904630 0.00146 0.00455 

Secondary event 

1409.998380 -0.00071 0.00455 



Table 2. Photometry for binary WTS 19g-4-02069 showing, from 
top to bottom, the WTS J-band photometry of WTS 19g-4-02069, 
and the INT i band data. Quoted magnitudes in the WFCAM sys- 
tem (column 3) can be converted to other photometric systems 
as described in Hodgkin et al. (2009). The J-band errors <7j, arc 
estimated using a noise model including Poisson noise, sky noise, 
readout noise and errors in the background estimation. (This ta- 
ble is published in full in the online journal and is shown partially 
here for guidance regarding its form and content). 

2.2 INT i'-band follow-up photometry 

We obtain follow-up photometric observations in the Sloan 
i-band on the 2.5m Isaac Newton Telescope (INT) on La 
Palma, using the Wide Field Camera (WFC), to refine our 
best-fitting light curve solution from the WFCAM J-band 
survey data. This imaging camera has a field of view of ap- 
proximately 34x34 arcmin 2 at prime focus, comprised of 
a mosaic of four 2 x4k pixel CCDs, with a resolution of 
~0.33 arcsec/pixel. The observations are part of a wider 
WTS follow-up campaign to confirm planetary transit can- 
didates, between July 18 and August 01 2010, leaving a few 
windows to observe binary eclipses. We use the WFC in fast 
readout mode (readout time 28 sec, for lxl binning) to ob- 
serve a full primary eclipse of WTS 19g-4-02069 on the night 
of July 25 2010 and a full secondary eclipse on the night of 
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Figure 1. The WFCAM J-band discovery light curve of binary WTS 19g-4-02069, before (left upper panel) and after (right upper panel) 
removal of the rotational signal, which reduces the total light curve rms by ~20%, and corrects the shape of the eclipses. The data are 
shown in relative magnitudes Am. The best-fit JKTEBOP models are overplotted with a solid red line. We also show zoom-ins around the 
primary and secondary eclipse to illustrate the effect of the correction. 



July 31 2010. We center the observations around the pre- 
dicted times of eclipse (based on the J-band data), and al- 
low ~50 min. of observations on either side of the predicted 
eclipses to acquire sufficient baseline. We obtain 82, 90 sec. 
exposures for the primary event and 52, 200 sec. exposures 



for the secondary, with out-of-eclipse rms of ~5.6mmag and 
~2.2mmag. Error-bars on the data are obtained by assuming 
a x 2 value of 1 for the out-of-eclipse parts using the models 
of Section 4.2. We estimate shot noise errors of ~3.3 mmag 
(90 s integrations) and ~2.2 mmag (200 s). 
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Figure 2. The follow-up INT i'-band light curves of the primary (left panel) and the secondary eclipse (right panel) of binary WTS 
19g-4-02069. Ovcrplottcd in red is the best-fit JKTEBOP model. 



We reduce the data using custom-built idl routines to 
perform the standard 2-D image processing (i.e. bias sub- 
traction and flat-field division). We remove low- level fringing 
by subtracting a scaled super sky-frame, which is obtained 
by median averaging ditthered exposures of a blank field un- 
der dark sky conditions. To generate the light curves, we use 
variable aperture photometry and circular apertures with 
the idl routine APER. We estimate the sky background us- 
ing a 3cr clipped median on 30x30 pixel boxes. A master ref- 
erence light curve is obtained from differential photometry 
on a set of ~10 bright, nearby, non-saturated, non-blended 
reference stars, selecting for each reference star the aperture 
that minimises the out-of-eclipse rms. Airmass dependence 
is removed by fitting a second order polynomial to the out- 
of-eclipse data. The INT i-band data are presented in Table 
2, and shown in Figure [2] 

2.3 WHT ISIS spectroscopy 

2.3.1 Low resolution reconnaissance spectrum 

We obtain low-resolution spectroscopy using the William 
Herschel Telescope (WHT), to confirm the M-dwarf nature 
of the binary system WTS 19g-4-02069 via measurement 
of the strengths of gravity sensitive atomic stellar absorp- 
tion lines, and to estimate the system effective temperature 
and chromospheric activity. Observations are carried out on 
the night of July 16, 2010 using the Intermediate disper- 
sion Spectrograph and Imaging System (ISIS). We use the 
R158R grating, which has a spectral resolution of R~1200 
(1.81 A/pixel) at 8500 A and wavelength range of ~6000- 
9000 A, for our 1.0" wide slit, chosen to match the typical 
seeing conditions of the night. We opt to use ISIS red arm 
only with the red sensitive RED+ array and not use the 
dichroic because it can cause systematics and loss of effi- 
ciency up to ~10%. A single 500 sec. spectrum is obtained 
using this setup at an airmass of ~ 1.35. 

We reduce the data using a combination of custom-built 
idl procedures and standard iraf routines. In idl we trim 
the spectrum, bias subtract and filter for cosmic rays, be- 
fore we divide our data by a median averaged flatfield, which 



we first correct for dispersion effects using a pixel-integrated 
sensitivity function. We use iraf's APALL routine to perform 
optimal 1-D spectral extraction. Wavelength calibration is 
obtained using arc spectra from the standard CuNe-CuAr 
lamps and flux-calibration and telluric line removal is at- 
tempted with an early-type standard star closely matched 
to the target in airmass. The final spectrum is shown in Fig- 
ure[3]as the black continuous line. In the upper panel of this 
figure, we indicate a few important molecular absorption 
bands, and several atomic lines. A best-fit model spectrum 
(solid red curve) and a best-matching observational template 
(solid green curve) are also overplotted on the lower panel 
of Figure [3] which are discussed in further detail in Section 
3.1. This spectrum is taken at 0.4 in binary phase, and al- 
though at this phase possible H Q emission of the secondary 
star may be blended with the primary we argue in Section 
5.3, using medium resolution H a spectral observations taken 
near quadrature, that such a contribution is likely very small 
and the shift of the H Q emission line is consistent with the 
primary star's radial velocity curve, indicating that the sec- 
ondary is very dim at these wavelengths and may not be a 
very active star. 



2.3.2 Intermediate resolution spectroscopy 

We obtain intermediate resolution spectra with ISIS on the 
WHT, in order to measure the binary's radial velocity curve. 
We use the red arm with the R1200R grating (spectral res- 
olution R~9300, 0.26 A/pixel) centred on 8500 A, giving a 
wavelength coverage of ~8100-8900 A. This wavelength re- 
gion (in particular the 8700-8850 A part) contains a number 
of relatively strong metallic absorption lines in M-dwarf at- 
mospheres. Using this setup we obtain 5, 60 min (3x20min) 
exposures on the nights of July 17 and July 18, 2010. Expo- 
sures of 90 min (3x30 min) and 85 min (1x1200+2x1800) 
are also acquired with the same grating, but centred on the 
H-alpha emission line (6563 A), on the night of July 28, 
2010. All observations are centered around the two quadra- 
ture points of the binary. Data reduction is performed us- 
ing standard iraf procedures for instrumental signature re- 
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Figure 3. Upper panel: The observed low-resolution spectrum of binary WTS 19g-4-02069 obtained with the WHT (solid black line), 
marking the main molecular and atomic features. The left most vertical dotted line marks H a in emission for WTS 19g-4-02069, suggesting 
active chromospheres. Lower panel: ISIS spectrum with overplotted the best-fit NexTGen model for T e yy=3200K (solid red line). The 
solid green line is an observed template spectrum for GJ1001, a single M4 star, archival data from the 6.5m Multiple Mirror Telescope 
(MMT, see Section 3.1). 
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moval (the CCDPROC package), with bias subtraction, flat 
fielding and correction for instrumental response. We then 
calibrate the observed wavelengths using CuNe+CuAr arc 
lamp spectra taken after each set of exposures. We flux cal- 
ibrate the data using observations of a photometric stan- 
dard. Our spectra have a relatively low signal-to-noise ratio 
of (S/N) = 15, meaning that because of the low luminosity of 
the secondary star in the optical, its spectral lines are not 
detected in this dataset. We use these observations to fur- 
ther constrain the radial velocity amplitude of the primary 
component as well as the orbital eccentricity in Section 3. 



2.4 GEMINI/GNIRS infrared spectroscopy 

We conduct observations with the GEMINI Near InfraRed 
Spectrograph (GNIRS) on the 8.1 m GEMINI-North Tele- 
scope, in queue-schedule mode in the H-band, to measure the 
secondary RVs. We move observations into the infrared, be- 
cause the cooler companion is brighter at these wavelengths 
and the estimated binary luminosity ratio in the infared ri- 
band is a factor ~3 higher than in the optical red spectra. 
We opt to use an intermediate resolution setup, rather than 
a high resolution setup, to maximise the spectral through- 
put for the secondary star, which results in a compromise 
on the velocity resolution. We use the long slit (49") and 
the long-red camera configuration, in combination with the 
110.5 lines/mm grating and a slit-width of 0.3", achieving a 
spectral resolving power of R~5900. Our set-up is centred on 
~1.555/xm, and has a wavelength range of 1.49-1.61/im. The 
corresponding velocity resolution is ~27.3 km/s/pixel. We 
observe our binary target using a standard ABBA on-source 
ditther pattern, nodding along the slit. In total, we obtain 
86, 240 sec. exposures on five seperate observing nights be- 
tween March and July, 2011. Of these 86 exposures, 79 have 
sufficient signal-to-noise for radial velocity work. The obser- 
vations are grouped as follows (corresponding to concurrent 
observations excluding low signal to noise data) : 6 x 240 s on 
the night of March 31, 2011, 12x240 s on May 31, 12x240 
s on June 1, 12x240 s on June 23 and finally 17 x 240s and 
21 x 240s on the night of June 17. All observing runs are 
centered around the two quadrature points of the binary 
(phases 0=0.25 and 0.75), where the relative RVs are ex- 
pected to be the largest. 

We reduce the data using the GNIRS sub-package (ver- 
sion vl.11.1) of the Gemini iraf package for spectral reduc- 
tion and extraction, which is available onlinqj, and which 
is adapted by us for optimal reduction of the current data 
set. Before running the iraf script we apply the python 
script CLEANIR.py to the raw data, to correct vertical strip- 
ing (repetative every 8 columns), horizontal banding, and 
quadrant offsets, which represents a significant source of 
additional background noise in ~70% of the 2-D spectra, 
the magnitude of which also varies between different data 
frames. In iraf, we first correct the data for read out noise, 
the detector offset (which is measured from a dark area 
of the data), and non-linear response using nsprepare. 
Using nsreduce we then remove the instrumental signa- 
ture by flatfielding and dark subtraction, followed by sky- 
subtraction by forming A— B and B — A pairs. On each sky- 



subtracted 2-D image we then measure the dispersion varia- 
tions using NSSDIST, and rectify the images using nstrans- 
form. With nscombine each A — B is then combined with 
the corresponding B — A by shifting the positive spectra 
on to each other, based on the header information. Because 
the raw data show frequent spikes from radioactive parti- 
cle hits, caused by decaying thorium on the lenses used in 
GNIRS, we median combine all sky-subtracted (A-B,B-A) 
pairs, which removes most of the hits. Our observations (es- 
pecially the part taken at high airmass) also suffer from 
highly variable sky-lines from OH sky-glow, which causes 
line residuals even after sky-subtraction. Optimal 1-D spec- 
tral extraction is then performed using nsextract, which is 
based on the standard iraf package APALL. Any remaining 
particle hits are identified by eye from the 1-D spectrum and 
clipped before further analysis. A summary of the spectral 
observations is shown in Table 3. 



3 SPECTROSCOPIC ANALYSIS 

3.1 Analysis of the low resolution ISIS spectrum 

We determine the absorption indices in the red part of the 
optical spectrum of Sodium (Nasi89) and Titanium Ox- 
ide (Ti07i4o) to constrain the luminosity class, and verify 
the main-sequence dwarf nature of WTS 19g-4-02069. We 
follow the procedure as outlined in Figure 11 of Slesnick, 
Carpenter & Hillenbrand (2006), to distinguish low surface 
gravity giants from high gravity dwarfs. We find indices 
Nagi89 ~0.9 and TiO7i40 ~1.9, which are consistent with 
a main-sequence dwarf of spectral type ~M3.5 (which is 
within the typical la uncertainties of the Reid, Hawley & 
Gizis 1995 spectral type relation for TiO). We use the metal- 
licity index Cmo /CaH (described by Reid et al. 1995, Lepine, 
Rich & Shara 2007 and Dhital et al. 2012), to find that the 
binary has solar metallicity within the uncertainties. Signif- 
icant H a emission is seen (we estimate an equivalent width 
EWfl =-6 A) which is probably caused by the magnetic 
activity of the binary, related to the chromospheres of the 
stars. The strong presence of the Na I doublet, the absense 
of deep infrared Calcium triplet absorption (8498, 8542 and 
8662 A), and lack of significant lithium absorption (6708 A), 
indicates that the primary is likely a mature M-dwarf, and 
not young and actively accreting nor a brown dwarf (e.g. 
Rebolo, Martin & Magazzu 1992). 

We follow the procedure outlined in Nefs et al. (2012) 
and Birkby et al. (2012), to determine the system effective 
temperature, which can be used to derive individual compo- 
nent temperatures Ti,2 when combined with the light curve 
parameters, by \ 2 fitting of a grid of NextGen atmosphere 
models (Allard et al. 1997) to the low resolution spectrum. 
This grid consists of models with AT e //=100K, and assumes 
constant log(g) = 5 (typical for field M-dwarfs) and solar 
metallicity, spanning the ~6000-9000 A region, which cor- 
responds to the data-range least affected by instrumental 
effects. The model spectra are scaled to match the contin- 
uum of the observed spectrum. In the fitting procedure, we 
mask the strong telluric Oxygen bandhead around 7600 A 
and the H a emission line. We use the formal errors as ob- 
tained from iraf to derive the x 2 , which we then optimise 



1 Available at http://www.gcmini.edu/sciops/instruments/gnirs/data-ffctmdilteriidHmid model. To derive an error, we 
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Table 3. A summary of the spectral observations obtained with WHT ISIS and Gemini GNIRS. For the R1200R ISIS setup, the 
spectrograph is centered on 8495 A, and around the H a emission line at central wavelength 6562 A. The columns RVi and RV2 indicate 
the radial velocity derived from the primary and secondary line shifts respectively, uncorrected for the systemic velocity of the binary 
system, but converted to the heliocentric system, as in Figure [4] 
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Figure 4. Radial velocity data for the primary (black filled symbols) and the secondary (red filled dots) stars of M-dwarf binary WTS 
19g-4-02069 using GNIRS. Black filled stars arc data from WHT ISIS, dots are from GEMINI GNIRS. The black and red solid curves 
are the best-fit sine functions. We obtain a binary systemic velocity of 29.1 (±0.5)km/s. The data in the figure are corrected to the 
heliocentric system. 
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scale to Xred = li yielding a final value of 140 K. Our best- 
fit model indicates T e //=3200 K, which is consistent with a 
~M3-4 spectral type, following the T e ff, spectral type re- 
lation from Baraffe & Chabrier (1996). Assuming that the 
primary star dominates the emission of the system in the op- 
tical, this spectral type roughly corresponds to that of the 
hotter component. The spectrum of GJ1001 (a ~M4 sin- 
gle nearby M-dwarQ), which was the best-fit observational 
template from a grid of M-dwarf spectra of various spec- 
tral subtypes observed with the Multiple Mirror Telescope 
(MMT), is overplotted on the lower panel of Figure [3] as a 
green continuous line. It is a better match to the ISIS data 
in the spectral regions around 6900 A than the NextGen 
model. 

We therefore identify WTS 19g-4-02069 as a genuine 
main sequence M-dwarf binary, with ~M3.5 and ~M5 com- 
ponents, and derive a best-fit effective binary temperature 
of 3200(±140)K. Its stars have nearly solar metallicity and 
the system is significantly active with an EWh q =-6 A. 



3.2 Radial velocities 

3.2.1 WET /ISIS 

We use the iraf routine FXCOR in conjunction with a grid 
of template synthetic stellar atmospheres of low-mass stars 
from the MARC^B spectral library (Gustafsson et al. 2008) , 
degraded to match the resolution of the observed data, to ob- 
tain radial velocities from the R~9300 ISIS spectra through 
1-D cross-correlation. In the cross-correlation procedure, we 
mask out the saturated near-infrared Ca II triplet lines at 
8498, 8542 and 8662 A. The spectral templates have a plane- 
parallel atmospheric geometry, an T e f / range of 2800-4000 K 
(in steps of 100 K), solar metallicity, log(p)=5 and 2.0km/s 
microturbulence. For the final RVs we use the template 
model that maximises the strenght of the cross-correlation, 
which is the cool T' e //=3200 K model. For the two H a obser- 
vations we simply fit single Gaussians to the emission line. 
The data is listed in Table 3 and plotted in Figure U as 
black filled stars. Although radial velocities measured from 
H-alpha lines can be biased depending on where the H-alpha 
emission region is located, we find that our two radial veloc- 
ity measurements are in good agreement with our red optical 
measurements and the line shape nearly Gaussian with no 
apparent major substructure. From the primary RVs we can 
already set reasonable upper and lower limits on the mass 
ratio of the system, because there is only a limited range of 
(Mi, Mi) that can yield the observed K\. We find an upper 
limit, assuming M\ =O.O8M0 (corresponding to the hydro- 
gen burning limit), of q <0.55, and a lower limit, assuming 
an MO-dwarf primary with Mi =0.65M Q , of q >0.25. This 
already indicates that WTS 19g-4-02069 is in the interesting 
regime of short period low-mass-ratio M-dwarf binaries. 

3.2.2 GEMINI/GNIRS 

We use the spectral region 1.55-1.6/im, which is the least 
contaminated by telluric water vapor absorption, to ex- 
tract RVs from the GNIRS data. With the IDL procedure 



Cross-correlation around phase 0.75 




Shift (km/s) 



Figure 5. The cross-correlation function for the summed spectral 
data for GNIRS Gemini around the 0.75 quadrature point using 
a template spectrum with T e f j-=3200K, corrected for the Solar 
motion and the systemic velocity of the binary. The red and the 
blue solid curves indicate the best-fit single Lorentzian profiles 
for the primary and secondary lines. 



cjcrosscorrelate and a grid of MARCS model spectra we ob- 
tain a cross-correlation function (CCF) for each spectrum, 
which we subsequently fit by a Lorentzian. To improve the 
contrast of the CCF, we first normalise the spectrum by a 
second order polynomial. We evaluate for each combination 
of template spectra for the binary components the corre- 
sponding strength of the CCF. We find that the T=3200 K 
model provides the highest signal for the primary of WTS 
19g-4-02069, whereas for the secondary a cooler template of 
T=3000 K maximizes its cross-correlation signal. Given the 
large orbital velocities, there is relatively little blending of 
correlation peaks. We report the detection of a clear sec- 
ondary component in the CCF around the <j> =0.75 quadra- 
ture point, shifted by ~4.5 spectral pixels. We show the CCF 
of the total added spectral data around <j> =0.75 in Fig- 
ure [5] (black curve) . The blue and red solid curves indicate 
best-fitting single Lorentzian profiles to the cross-correlation 
peaks of the secondary and primary binary components re- 
spectively. The data around the 0.25 quadrature suffer from 
higher noise due to highly variable sky and detector pattern 
noise, making it difficult to resolve the secondary lines. Only 
the datapoint at phase 0.28 has sufficient quality for radial 
velocity work on the secondary. 

We fit the RVs as function of binary orbital phase <f> with 
a simple sine curve using idl's mpfitfun function, assuming 
a circular orbit. We first fit the primary RVs from ISIS and 
GNIRS, allowing only K\ and 7 (the systemic velocity of the 
binary) to vary, fixing the phase using the well-determined 
orbital period from the light curve. For the secondary star 
we fix 7 to the value derived from the primary and fit for 
K2 (note that we obtain consistent results when we treat 
7 as a free parameter). To derive the RV errors we scale 
the errors from the Lorentzian/Gaussian CCF fit, such that 
the reduced \ 2 °f the best-fitting RV model is unity. We 
show our results in Figure [4] where the black solid curve 



indicates the best-fitting sine function to the primary RVs, 

2 Available at |http://spider.ipac.caltech.edu/staff/davy/ARCHIVE/ ind»tci;h1tM red solid curve the fit to the secondary RVs. We 

3 Available at |http: / /marcs. ast ro.uu.se71 find radial velocity amplitudes of A'i=29.4(±0.5) km/s and 
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Parameter 


i-band 


J-band 


J-band 

(spot-corrected) 


light curve analysis 








P (days) 


2.44178(±0.00003) 






T (MJD) 


2454319.83270(±0.00002) 






(Ri +R 2 )/a 


0.103(±0.002) 


0.104(±0.002) 


0.098(±0.004) 


R2/R1 


0.338(±0.003) 


0.343(±0.003) 


0.339(±0.005) 


J 


0.40(±0.01) 


0.65(±0.01) 


0.64(±0.02) 




87.7(±0.1) 


87.7(±0.2) 


88.2(±0.3) 


Ri/a 


0.077(±0.001) 


0.077(±0.001) 


0.073(±0.003) 


R 2 /a 


0.026(±0.001) 


0.026(±0.001) 


0.025(±0.003) 


L2/L1 


0.042(±0.004) 


0.074(±0.004) 


0.072(±0.007) 


T 2 /Ti 


0.779(±0.008) 


0.892(±0.007) 


0.89(±0.01) 


V 2 


1.093 


4.057 


1.004 


Estimated temperatures 










3200(±140) 






Ti(K) 


3300(±140) 






Ti(K) 


2950(±140) 






Radial velocity analysis 








R"l(km/s) 


29.4(±0.5) 






if 2 (km/s) 


109.0(±1.6) 






7(km/s) 


29.1(±0.5) 






Q 


0.27(±0.02) 






a(il Q ) 


6.7(±0.2) 






Derived masses and radii 








A/i (M ) 


0.53(±0.02) 


0.53(±0.02) 


0.53(±0.02) 




0.51(±0.01) 


0.52(±0.01) 


0.49(±0.02) 


M 2 (M ) 


0.143(±0.006) 


0.143(±0.006) 


0.143(±0.006) 


R 2 (R ) 


0.174(±0.006) 


0.177(±0.006) 


0.167(±0.009) 



Table 4. Best-fit parameters and derived quantities for the M-dwarf binary system WTS 19g-4-02069. The boldface fonts indicate the 
adopted masses and radii derived from the INT i-band light curve parameters (see the discussion in Section 4.3). 



#2 = 109. 0(±1. 6) km/s, for a systemic velocity 7=29.1(±0.5) 
km/s, indicating a binary mass-ratio g=0.27(±0.02). 



4 LIGHT CURVE ANALYSIS 

4.1 J-band photometry baseline variation analysis 

There is significant out-of-eclipse scatter in the J-band light 
curve (~2% peak-to-peak), which is not in phase with the 
binary orbit. To investigate the possible periodicity of this 
variation, we clip the primary and secondary eclipses from 
the light curve, and perform a frequency analysis using the 
idl implementation fasper.pro of the Lomb-Scargie peri- 
odogram (Lomb 1976; Scargle 1982; Press & Rybicki 1989). 
The left panel of Figure [6] shows this frequency spectrum 
for binary WTS 19g-4-02069. Using the idl routine, we de- 
termine a false-alarm probability (FAP) to filter out peaks 
that are likely caused by spurious detections on light curve 
systematics (horizontal dashed line). Significant power is ap- 
parent around ~2.44 d binary period, yet the actual peak is 
~0.14 d away at ~2.56 d. We indicate various integer and 
half-integer aliases of the binary orbital period in Figure 
[5] using vertical dotted lines. From previous binary stud- 
ies, similar deviations are seen in systems that are either 
young and not fully synchronised (e.g. the pre-main sequence 
eclipsing binary Paranengo 1802; Cargile et al. 2008) or not 
fully circularised due to their relatively long orbital period 
(e.g. Irwin et al. 2011). We discuss possible causes for this 
apparent discrepancy in Section 5. 

When folded onto the non-synchronous 2.56 d period 
determined from the Lomb-Scargie analysis, we find a clear 



nearly-sinusoidal modulation of the data, which we at- 
tribute to star spots on a rotating stellar surface. We also 
perform an independent check of these results using the 
idl implementation epfold of the Analysis of Variance 
(AoV;Schwarzenberg-Czerny 1989) algorithm and plot the 
results in the right panel of Figure [6] We confirm a best-fit 
period which is very close to the ~2.563 d period suggested 
by the Lomb-Scargie algorithm. No significant signal is ob- 
tained at the orbital period. The second highest peak in the 
right hand panel is at ~2.432 d, but at relatively low signif- 
icance. In Figure [7] we show the clipped light curve folded 
on the best-fit periodicity for spot modulation. We attempt 
to correct the J-band light curve using a single sine with 
8.1 mmag amplitude, which appears to minimise the out-of- 
eclipse rms. In the upper left panel of Figure [1] we show the 
phase-folded but uncorrected light curve (black filled dots), 
while in the upper right panel panel we show the J-band 
data, corrected for the ~8.1mmag rotation signal. We do 
not claim that this method is the best method for remov- 
ing rotational signal from a light curve, because if spots are 
occulted during the eclipse the light variation may strongly 
deviate from the sinusoid. However, we find that with our 
correction the out-of-eclipse rms decreases from ~10.6 to 9.1 
mmag (~15% reduction), and in eclipse it decreases from 
~11.2 to 8.9 mmag (~20% reduction). This means that the 
rms in and out of eclipse are approximately the same after 
the correction (in fact the difference between the rms values 
in and out of eclipse decreases from to ~2% after cor- 
rection). Zoom-ins around the eclipses of the corrected and 
uncorrected data are shown in Figure [T] 
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Figure 6. Left panel: Lomb-Scargle diagram of the dominant frequencies in the WFCAM J-band light curve after removing the eclipses. 
A significant peak is at ~2.563 days. We indicate several period aliases with the red dashed lines. The red horizontal line indicates the 
99% confidence level for peak rejection. Right panel: Analysis of Variance results for binary WTS 19g-4-02069. We show the frequency 
spectrum between 2.1 and 2.7 d and note that we confirm a strong signal at the period determined by the Lomb-Scargle method, and 
no significant signal at the orbital period P or h ~2.44 d. The second most significant peak is at ~2.432 d, but with low significance. 



4.2 jktebop light curve fits 

4-2.1 J-band photometry 

We use the binary light curve modelling program jKTEBOfQ 
(Southworth et al. 2004), which is based on the Eclipsing 
Binaries Orbit Program (ebop; Popper &Etzel 1981; Etzel 
1981) for the fitting of the uncorrected and corrected J-band 
data. For binary WTS 19g-4-02069, with orbital period 2.44 
d and mass ratio 0.27, Equation 6 of Morris (1985) predicts 
ellipsoidal light variations of just ~0.4 mmag in the J-band, 
indicating that the binary stars are likely only very slightly 
deformed by mutual tidal interactions. This, together with 
the low derived values of the stellar oblateness in our subse- 
quent fitting with jktebop, justifies the application of this 
model (which is only suitable for detached systems) to the 
WTS 19g-4-02069 binary system (see also Popper & Etzel 
1981). 

For the light curve modelling we allow the following six 
parameters to vary: i) the binary orbital period Po, ii) the 
mid-eclipse epoch To of the primary eclipse, iii) the sum 
of the stellar radii in units of the binary semi-major axis, 
(Ri + R,2)/a, iv) the ratio of the radii k = R2/R1, v) the or- 
bital inclination i and vi) the ratio of central surface bright- 
ness J = J2/J1. We use as input the initial estimate of the 
binary ephemeris obtained from the automated OCCFIT BLS 
algorithm. 

We keep the mass-ratio of the system, q, fixed at the 
value determined by the spectroscopic analysis. We do not 
fit the reflection coefficients, but calculate them from the 
system geometry. We also assume a gravity darkening coef- 
ficient, which is fixed at a value typical for stars with con- 
vective envelopes (/3=0.32; Lucy 1967). We find that treat- 
ing additional third light as an extra fitting parameter does 
not significantly improve the quality of our fit, so we fix 
it to zero. We adopt linear limb-darkening coefficients from 

4 Available at |http: / /www. a stro, keele,ac,u k/~jkt/| 



t-band J-band 
Linear Primary 0.69 0.41 
Secondary 0.78 0.46 
Quadr. Primary [0.22;0.58] 
Secondary [0.37;0.49] 

Table 5. Limb-darkening coefficients used as input to the EBOP 
models for the WFCAM J-band and the INT i-band light curves. 
For the optical red data we indicate both linear and quadratic 
coefficients. 



Claret et al. (2000) in the J-band (see Table [SJ , which are 
calculated from PHOENIX model atmospheres (Allard et 
al. 1997), for a surface gravity log(g) = 5, solar metallicity 
and 2 km/s micro-turbulence and stellar effective tempera- 
tures 21,2, such as derived in Section 4.2.2. We do not fit for 
the limb-darkening coefficients, because the S/N is too poor, 
and keep Ti^ (see Section 4.2) fixed. The orbital eccentric- 
ity e and the argument of periapse ui are also kept fixed to 
values consistent with a circular orbit, because our initial 
runs indicate that the data are firmly consistent with such 
an orbit (|e*cos(cj)| < 0.000079). This is expected from tidal 
dissipation theory given the relatively short circularisation 
timescale ~200 Myr (Zahn 1977). 

To assess the la parameter uncertainties we use the 
Monte Carlo routine from jktebop (Southworth et al. 
2005). In this procedure, Gaussian random noise is repeat- 
edly (10000 times) added to the model light curve before a 
new fit is made to the data, which yields a distribution for 
each parameter. With jktebop we also perform a prayer- 
bead error analysis, which can be useful in the presence of 
correlated noise (Southworth 2008), and find that parameter 
values derived from both methods are consistent within the 
la uncertainties, and that the uncertainties are not signifi- 
cantly different from the prayer-bead analysis. We therefore 
adopt the MCMC method from hereon. Numerical results 
of the light curve fitting are given in Table 4. 

The derived parameters for the uncorrected and the 
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Figure 7. The clipped WFCAM J-band light curve folded on the 
~2.563 day frequency peak. A significant rotational modulation 
can be seen. The solid red curve is the best fit sinusoid to the 
data with amplitude ~8.1mmag. 



spot-corrected data are inconsistent within the quoted la 
uncertainties, (Ri + 7l 2 )/a, R2/R1, and J are smaller after 
correction, whereas i is higher. In Figure [T] it is apparent 
that the correction removes outlying data points in eclipse, 
which bias the measured depths of the eclipses and the dura- 
tion of ingress/egress. Note that the errors of the corrected 
light curve data points are scaled to obtain X^ ed —1. The 
derived ratio of radii (i?2/i?i=0.339±0.005) in the data in- 
dicates that the companion is significantly smaller than the 
primary. An impact parameter b = (a/Ri) *cos(i)=0A3 sug- 
gests that the eclipses are full and the system is non-grazing, 
which lifts the degeneracy between R2 / Ri and i. The ratio of 
secondary to primary luminosity is Z/2/£i=0.072(±0.007), 
which indicates that only ~7% of the system light in the in- 
frared is due to the secondary star. From L2/L1 and R2/R1 
we derive, using Stefan-Bolzmann's law, a wavelength spe- 
cific temperature ratio T2/Ti=0.89(±0.01). Assuming that 
all of the light in the system comes from the primary, and 
the stars radiate as blackbodies, this indicates T2 ~2850K 
for the secondary, which would be consistent with a M5 type 
star according to the T e //-spectral type relation presented 
in Stephens et al. (2009), M5 according to the 1 Gyr model 
from Baraffe & Chabrier (1996), and M6 according to Reyle 
et al. (2011). The mass and radius of the secondary, such as 
derived in Section 4.3, indicate a surface gravity log(g)~5.1 
which is consistent with a main-sequence M-dwarf. 



4-2.2 i'-band photometry 

The lack of out-of-eclipse baseline for the optical data means 
that we can not accurately determine the amplitude of any 
spot modulation, so we opt to fit the light curve with- 
out making any spot corrections. The INT eclipses have 
a flat bottom, which confirms that the secondary star is 
fully superimposed on the primary during eclipse. To model 
with jktebop, we adopt the linear coefficients from Claret 
et al. (2004) for the Sloan i'-band. The derived parame- 
ters (Ri + -f?2)/a and R2/R1 are consistent with the spot- 
corrected infrared results, within the quoted la uncertainties 



of the J-band data. The orbital inclination i is slightly lower 
in the optical, but only by ~1.4er. The optical data reveal 
a surface brightness ratio J=0.40(±0.01) and a wavelength 
specific luminosity ratio L 2 /Li=0.042(±0.004). This result 
shows that the secondary is significantly dimmer at optical 
wavelengths. Also, a ratio T 2 /Ti=0.779(±0.008) is derived 
in the optical (~11% lower than in the infrared). To derive 
component temperatures Ti.,2 we use MARCS model spec- 
tra, convolved to the i' and J bandpasses, to derive model 
surface brightness ratios J\, which are compared to the op- 
tical and infrared observations. We reproduce the observa- 
tions, within the estimated la uncertainties, for Ti=3300 
K and T2=2950 K , with an estimated uncertainty of 140 
K. We adopt these values of Ti and T2 in our subsequent 
discussion. 



4.3 Stellar Masses and Radii 

We derive the component masses and radii from the com- 
bined RV analysis (incorporating both the ISIS and the 
GNIRS measurements), the i-band light curve fitting pa- 
rameters, and Keplers law according to the following two 
equations (see e.g. Hilditch 2001): 



Mi = 



2iYGsin(i) 3 



Ri 



Xia ■ 



P 2 GM!(1 + K 1 /K 2 ) 



47T 2 



1/3 



(1) 

(2) 



where x\ denotes the best-fit scaled binary orbital separa- 
tion, Ri/a, from jktebop. We propagate the errors from 
the light curve and radial velocity analysis. The main mo- 
tivation for using the i-band results is that for the infrared 
data, many spot cycles are folded into the light curve, and 
although our single sine correction removes part of the scat- 
ter, tracing the stellar activity cycle over such a long ob- 
servational baseline (~5 years) is difficult because the spot 
configuration may have evolved significantly. Moreover, the 
J-band data has significant gaps in observing dates and gen- 
erally only few observations per night, making it difficult to 
accurately model the precise behaviour of the spots. Be- 
cause the optical data for the primary and the secondary 
event are obtained within 6 d of observations, they likely 
trace the same spot configuration. Moreover, both the pho- 
tometric quality and the number of in-transit datapoints 
of the optical data can rival the J-band photometry. One 
problem in the current work is the limited coverage of 
the parts out of eclipse. Work by Goulding et al. (2012) 
shows that the light curve amplitude of a spot in the J- 
band is generally 55% of that in the I-band, indicating 
that ~3% peak-to-peak variations would be expected for 
INT observations of the full binary orbit. This could intro- 
duce an additional error to our light curve fitting results. 
Future work should adress this issue by obtaining multi- 
band photometric observations of concurrent binary eclipses 
and sufficient baseline, to catch a single spot cycle. With 
2i=0.077, A'i=29.4km/s, and if 2 =109.0 km/s, we find radii 
of i?i=O.51(±O.Ol)i? , i?2=O.174(±O.OO6)7? , and masses 
Mi=0.53(±0.02)M Q , M 2 =0.143(±0.006)M Q . We use equa- 
tions similar to Equations 1 and 2 to derive the mass and 
radius of the secondary. This translates to 2.1-3.2% errors 
on the radii and ~4% errors on the masses, although we 
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caution that these errors do not include possible uncertain- 
ties from star spots. Note that these masses and radii are 
consistent with main-sequence model predictions (see Sec- 
tion 5.2), rather than pre- main sequence, providing further 
support for the mature nature of the system. 



5 DISCUSSION 

In this paper we have presented the discovery of a highly 
unequal-mass eclipsing M-dwarf stellar binary, using the 
high-precision infrared light curves of the WFCAM Transit 
Survey, and follow-up characterisation with optical photom- 
etry and optical and IR spectroscopy on 2.5-8 m class tele- 
scopes. With two components straddling the fully convective 
boundary, and with shared ages and metallicities, our binary 
provides a rare and more stringent comparison to structure 
model predictions of fundamental M-dwarf properties over 
a wide span of stellar masses. The cool M5V secondary of 
the binary is in an important mass-regime for studies of 
exoplanets down to the Earth-size regime. In this section 
we will discuss our binary in the context of current theo- 
ries for low-mass binary formation, which predict such close 
unequal systems to be rare and in the context of low-mass 
stellar structure models. 

5.1 The mass-ratio distribution 

The distributions of binary orbital separation and mass- 
ratio, as function of primary mass, provide important con- 
straints on star formation simulations (e.g. Burgasser et al. 
2007, Bate et al. 2012, Clarke 2012). These simulations sug- 
gest that accreting gas with high angular momentum and 
dynamical interactions tend to drive up the mass-ratio of 
close binary systems towards unity. Dynamical interactions 
are frequent in the binary birth environment, and inter- 
actions with more massive stars will generally bias binary 
primaries towards higher masses. This indicates that it is 
unlikely that close, unequal binary systems with low-mass 
primaries can be maintained for very long. For example, a 
0.5+0. IMq M-dwarf binary that is formed in a cluster with 
stellar density n, ~2000pc~ 3 has a life-expectancy against 
disruption from solar type stars of order ~10 Myr (Goodwin 

6 Withworth 2007). Furthermore, low-mass binary systems 
and unequal-mass systems are more easily perturbed due to 
their lower binding energy. 

In Figure [8] we show the mass-ratio distribution of close 
binary systems (P or b <10 d) with M-dwarf primaries. This 
figure is compiled from Table 6, which holds the currently 
available sample of M-dwarfs discovered as eclipsing and 
non-eclipsing, double-lined, spectroscopic binaries. This is 
an updated version of Figure 9 in Wisniewski et al. (2012), 
who include only a few M-dwarf systems in the eclipsing bi- 
nary period range. The left panel of our Figurc[5]shows mass- 
ratio as function of binary orbital period for 55 sources, of 
which 20 are spectroscopic binaries, whereas the right panel 
shows a histogram of the mass-ratio distribution in bins of 
0.05. We separate the M-dwarfs depending on whether their 
primaries are more or less massive than the fully convec- 
tive boundary at ~O.35A/0 (black and red filled squares). 
With a mass-ratio of 0.27, it is clear that WTS 19g-4-02069 
occupies an interesting position in these diagrams, because 



over 80% of the stellar binaries have q>0.8. There may be 
an observational bias towards more equal-mass M-dwarf bi- 
naries because of the steep relation of mass and luminosity 
for M-dwarfs, causing the spectral lines of low-luminosity 
companions to remain unresolved with optical spectroscopy. 
Also, for lower q systems, lower radial velocity shifts of the 
primary lines are expected, making them more difficult to 
detect a priori, especially for binaries with late type M-dwarf 
primaries, which are intrinsically faint. In order to better 
understand this bias it is important to move future spectro- 
scopic and photometric observations into the infrared, be- 
cause of the improved binary luminosity ratio there. Also 
any single- lined short-period eclipsing systems should be re- 
ported, because for such systems an upper limit can be ob- 
tained on the mass-ratio given that the orbital inclination 
can be directly constrained from the light curve. Yet, note 
that in the Birkby et al. (2012) candidate MEB list, more 
particular the sample with best-fit binary temperature less 
than 4000 K, we do not identify any object with a ratio of 
secondary to primary eclipse depth significantly lower than 
that of binary WTS 19g-4-02069, even though we have the 
photometric precision to detect such low-mass M-dwarf com- 
panions through their eclipses. This may indicate that such 
systems are in fact intrinsically rare. As we will discuss in 
the next paragraph such considerations are important con- 
straints on the physics of low-mass binary formation and gas 
accretion mechanisms. 

Three binaries are observed in the range q=[0.4-0.5], 
but none at lower q. Of these three, two are very short- 
period (~0.4 d) and young, 20 Myr T-Tauri stars (NSVS- 
06507557; Cakirh & Ibanoglu 2010) and 150 Myr young clus- 
ter members (2MASSJ04463285; Hebb et al. 2006). Theory 
shows that significant dynamical processing can occur prior 
to the main-sequence, indicating that low mass-ratio sys- 
tems should be more abundant while young. For example, a 
significant difference in the binary fraction between young 
clusters and field solar type stars has been observed (e.g. 
Duchene et al. 2007). 

A possible explanation for the existence of WTS 19g-4- 
02069 is that the physics of gas accretion onto (close) binary 
systems is different than suggested by smoothed particle hy- 
drodynamic (SPH) and ballistic particle simulations (from 
e.g. Artymowicz 1983; Bate & Bonnell 1997; Bate 2002). 
Two-dimensional warm grid-based simulations from Ochi et 
al. (2005) and de Val-Borro et al. (2011) confirm that gas 
preferentially enters the secondary Roche lobe, but flows 
around the secondary and is then channeled onto the pri- 
mary star, which grows a more massive accretion disk, which 
means that q can decrease during the accretion phase. Very 
recent numerical work from Zhao & Li (2012) suggests that 
adding magnetic fields to binary formation simulations can 
have a similar effect. Such a magnetic field could apply a 
brake on the material that flows onto the binary, decreas- 
ing its angular momentum and significantly shrinking the 
protobinary separation, meaning that q could be low for 
short binary orbits. If these simulations are correct, low- 
q pairs are expected to be abundant around a wide range 
of primary masses on the main-sequence, which would be 
consistent with our observations of WTS 19g-4-02069. How- 
ever, disrupting third-body interactions can still remove low- 
q components from binaries. Another explanation for the low 
mass-ratio of WTS 19g-4-02069 could be that the binary was 
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Figure 8. The mass-ratio distribution of short period (P <10 d) M-dwarfs reported in literature. Left panel: binary orbital period (in 
d) versus mass-ratio. Black filled squares indicate eclipsing M-dwarf binaries with Mi > O.35M0, red squares indicate eclipsing binaries 
with Mi < O.35M0, whereas open diamonds show M-dwarf binaries reported in radial velocity surveys. The large blue cross is our 
measured mass-ratio for M-dwarf binary WTS 19g-4-02069. The three black dashed lines show the respective mass-ratios for M-dwarf 
binary systems with a secondary at the hydrogen burning limit, assuming primaries of 0.6A/q (M0 spectral type; Baraffe & Chabrier 
1996), 0.2Mq (M3), and O.IIMq (M5). These lines indicate the lower limits in q to which M-dwarf primaries with stellar secondaries 
are confined. Right panel: histogram of mass-ratio. 



isolated from the birth environment early on, for example 
through ejection from the cloud due to binary-single star 
or binary-binary interactions, or because the natal cluster 
became unbound due to rapid gas removal. 

5.2 The mass-radius relation for M-dwarfs 

Highly unequal-mass M-dwarf binaries provide important 
test cases of low-mass stellar evolution theory because they 
cover a large range of M-dwarf masses, which encompass 
significant changes in stellar atmospheric structure. Further- 
more, since the binary components have the same age and 
metallicity, stellar evolution models can be tested using two 
less free parameters. Also, Earth-like planets in the habit- 
able zones around late-type M-dwarfs are prime targets for 
new transit surveys and Extremely Large Telescopes (ELTs), 
and require accurate calibration of their small host stars. 
For example, for the Neptune-sized planet orbiting the M5V 
host star GJ1214b (Charbonneau et al. 2009, de Mooij et al. 
2012), a 15% uncertainty in stellar radius could translate to 
the difference between an ocean planet and a gaseous H- or 
He-type atmosphere. 

Theory predicts that stars in the fully convective mass 
regime (<O.35M0) respond differently to rapid rotation and 
strong magnetic fields (which can reduce convective effi- 
ciency) than their partially-convective cousins (e.g. Chabrier 
et al. 2007). Most importantly, because in fully convective 
atmospheres heat flow is nearly adiabatic, model stellar radii 
and temperatures are expected to better match the observa- 
tions. It is argued that strong magnetic fields can shift the 
fully convective mass boundary to masses as low as O.IMq 
(Mullan & MacDonald 2001). Such predictions can be tested 
with a sample of mass-radius-temperature measurements. 
However, in the stellar-mass regime <O.2M0 there is a lack 
of model-independent data. 

In Figures [9] and [10] we present the current sample 



of mass-radius-temperature measurements for M-dwarfs. In 
these figures, filled black and red dots indicate dynami- 
cally derived model-independent masses and radii for the 
primaries and secondaries of eclipsing double-lined M+M- 
dwarf binary systems. Green dots represent other measure- 
ments of M-dwarfs which mostly rely on model-dependent 
assumptions and/or external constraints on the system. 
These systems include M-dwarfs as secondaries of single line 
binary systems, tidally interacting white dwarf- M-dwarf sys- 
tems and single star M-dwarf measurements which often as- 
sume an empirical mass-luminosity relation. In the inset of 
Figure [9] we show the available measurements in the regime 
<O.2A/0, further detailing their origin. There is significant 
scatter in these data. Four of the data-points are from M- 
dwarfs orbiting F- or G-type stars, i.e. single line systems, 
and rely on model-dependent constraints on the properties 
of the primary and/or assume spin-orbit alignment (Pont 
et al. 2005, 2006; Beatty 2007). Three of the M-dwarf sys- 
tems have either a white dwarf primary, which may have had 
a phase of common envelope evolution or significant mass- 
transfer and are likely to be tidally interacting, or B-type 
subdwarf (Parsons et al. 2012;Pyrzas 2012). Three other 
systems are from interferometric data with directly mea- 
sured radii, but with estimated masses from a model mass- 
luminosity relation (e.g. Demory et al. 2009). Finally two 
data-points are pre-main sequence M-dwarf binaries which 
have secondaries in the <O.2M0 regime (JW380, Irwin et 
al. 2007; 2MASSJ04463285, Hebb et al. 2006). The cool 
companion of WTS 19g-4-02069A is in a special position 
as it provides a model-independent anchoring point of mass 
and radius for fully convective low-mass main-sequence M- 
dwarfs. 

The companion is O.OO67R.0 (4.0%) larger in radius than 
the lGyr solar metallicity Baraffe model, whereas the pri- 
mary is larger by O.O25R.0 (5.0%). The radius of the compan- 
ion can be constrained to O.OO56i?0 (3.2%) at the la level, 
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Figure 9. The mass-radius diagram for low-mass stars with mass less than 0.7 Mq. The vertical dashed line at 0.075Mq indicates 
the hydrogen burning limit, whereas the line at <-->0.35Mq represents the proposed transition to fully convective stellar atmospheres 
(Chabrier & Baraffe 1997). Black filled dots are the primary stars in double-lined eclipsing M+M-dwarf binaries, whereas red filled dots 
are the secondaries. The small filled green dots indicate M-dwarfs with masses and radii determined through other methods, which we 
further examine in the inset for very low mass M-dwarfs (M<0.2Mq). The purple, blue, and yellow triangles indicate measurements of 
M-dwarfs in eclipsing systems with higher mass G or F type primaries, white dwarf (WD) primaries, or subdwarf B-typc (SdB) primaries 
respectively. The green filled triangles are single star measurements using radius measurements from interfcrometry. 



its mass is uncertain by O.OO56M0 (4.0%), indicating that it 
is currently outside the <3% range advocated by e.g. Tor- 
res (2012) as a stringent constraint on models. We illustrate 
this by showing three Baraffe (1998) models for: i) 1 Gyr 
solar metallicity (solid black curve); ii) 5 Gyr solar metal- 
licity (dotted curve) and hi) 1 Gyr metal poor, [M/H]=-0.5 
(dashed curve) in Figure [9] A 3% accuracy range on the 
companion mass could be reached if the uncertainty on K2 
is pushed down to the ~1 km/s level. It is interesting to 
note that the effective temperature for the fully-convective 
secondary is consistent with the 1 Gyr model within the 140 
K measurement error, but the partially-convective primary 
has a T e ff lower by ~350 K. Although there is significant 
scatter in the mass-T^// diagram around 0.5Mq, potentially 
due to the inhomogeneous set of methods by which temper- 
ature is determined, such a difference could be explained by 
magnetic inhibition theory (e.g. Chabrier et al. 2007) if the 
primary is an active star and has a significant magnetic field, 
whereas the secondary may not. 

Five of the seven literature eclipsing binaries with q<0.6 



are fast rotators (P <1 d), and three of these seven are still 
on pre-main sequence tracks. The primary of main-sequence 
M-dwarf binary V405 And (Vida et al. 2009; q ~0.429, 
P ~0.465 d) is inflated with respect to the Baraffe model 
by 73%, whereas the secondary is not significantly inflated. 
This system is much more active than WTS 19g-4-02069, 
with frequent flaring events. In NSVS-65550671 (Dimitrov 
& Kjurkcieva 2010, q ~0.510, P ~0.193 d), both compo- 
nents are inflated (~17% for the primary, ~12% for the sec- 
ondary). Hebb et al. (2006) present a q ~0.404, P ~0.465 
d system in the young open cluster NGC 1647, which has 
an inflated O.47M0 primary and a non-inflated O.19M0 sec- 
ondary. According to Figure 1 from Irwin 2007, 0.5Mq M- 
dwarfs reach the main sequence by ~150 Myr, and by ~300 
Myr for O.2M0, indicating that one or both components of 
the Hebb et al. binary may be pre-main-sequence. In the pre- 
main-sequence binary NSVS-06507557 (Cakirh & Ibanoglu 
2010; q ~0.425, P ~0.520 d) the contrary is true when com- 
paring with the models: the secondary is inflated by 35%, 
whereas the primary is consistent with the models. One ex- 
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Figure 10. Current census of mass, effective temperature measurements for stars with mass less than 0.7Mq. Similar to Figure[9] filled 
dots indicate primary (black) and secondary (red) components of double-lined eclipsing binaries. The blue crosses indicate the position 
of the components of binary WTS 19g-4-02069. The solid black, dotted and dashed curves are Baraffc model predictions for different 
ages and metallicities. The data for this figure are tabulated in the Appendix. 



planation is that in such a young system, the size contrac- 
tion through gravitational collapse of the secondary may not 
have progressed as much as that of the higher mass primary, 
because lower mass stars reach the Zero Age Main Sequence 
(ZAMS) at later times. Binary WTS 19g-4-02069, for which 
we find no clear indications of a young age, has a >5 times 
lower rotation rate than these four systems, suggesting that 
inhibition by stellar rotation is the main reason for the ob- 
served difference in radius inflation. 

5.3 Non-synchronous rotation? 

In Section 4.1 we determine that the WFCAM J-band data 
exhibit a ~2% peak-to-peak near-sinusoidal periodicity. The 
variability is not in phase with the binary eclipses, there- 
fore we attribute it to a rotational modulation in brightness 
caused by star spots. The H n profile, as observed with ISIS 
(Section 3.2.1), suggests that the primary star is the main 
contributor to both the emission of this line and the J-band 
variability, because the observed radial velocity shift is con- 
sistent with the primary. Furthermore, the secondary star 
contributes only ~7% of the total J-band emission, there- 
fore requiring spot modulations of order ~0. 0162/0. 07~0. 23 



mag, to account for the variability, which would require a 
very high spot coverage on the secondary and a magneti- 
cally quiet primary. The best-fit period of ~2.563 d, ~0.14 
d (~5%) longer than the binary orbital period of ~2.44 d, 
would suggest that the primary star rotates at a subsyn- 
chronous rate. This finding is in contrast with current pre- 
dictions from tidal theory (e.g. Zahn 1977), that suggest 
orbital synchronisation on timescales of the order 10 4 ~ 6 yr 
and sitpereynchronous rotation for young stars. We discuss 
two possible scenarios that could explain the observations. 

In the first scenario, high latitude spots have a longer 
rotation period than spots near the M-dwarf equator due 
to significant differential rotation over the stellar surface. 
This causes the observed rotation rate Q bs = 2n /Pots to 
be lower than the true equatorial rate fl eq . On the Sun for 
example, a spot at 60° latitude has a rotation period ~25% 
longer than on the equator, and the intermediate latitude 
6 is given by the equation Q(9) = Q. eq — dQ.sin 2 {6), where 
df2 is the difference in rotation between the equator and 
the poles. Donati et al. (2008) present a spectropolarimet- 
ric survey of a small number of single early-type M-dwarfs 
(ranging from M0 to M3), and detect significant differen- 
tial rotation in four objects of their sample. For example, 
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OT Ser (Ml. 5 dwarf; ~3.38 d - dQ./Q, £q ~0.06) and DS 
Leo (MO dwarf; ~14 d - dQ./Q. eq ~0.16) show surface ro- 
tation variations which are consistent with, although some- 
what higher than our observations. Barnes et al. (2005) ar- 
gue that for low-mass stars the surface differential rotation 
vanishes with increasing convective depth, such that fully 
convective stars rotate mostly as solid bodies. This is sup- 
ported by observations from e.g. Donati et al. (2006a) and 
Morin et al. (2008a,b) who show little to no differential ro- 
tation in late-type M-dwarfs. We therefore argue that the 
variability of WTS 19g-4-02069 could be caused by a large 
high-contrast stable cool spot complex located near to the 
rotation pole(s) of the 0.5Mq primary, where the rotational 
shear is relatively low. Similar large and stable (nearly) po- 
lar spot patterns with (quasi-)sinusoidal light curves are ob- 
served on the young ~M2V star AU Mic (see e.g. Rodono et 
al. 1986). A dichotomy in magnetic field geometry is pointed 
out by Morin et al. (2010) for stars above and below 0.5Mq, 
suggesting that the dynamo mechanism for 0.5Mq stars is 
rather like that of solar-type stars, where the Coriolis force 
induced by fast stellar rotation (WTS 19g-4-02069 rotates 
at ~10x the solar rotation value) tends to drive up star 
spots to polar latitudes. The photometric signal in Figure [7] 
has a relatively high scatter compared to the sine amplitude 
(~10.4 mmag scatter for a ~16.2 mmag signal), which may 
suggest that the size of the region could have fluctuated over 
the timescale (~5 years) of our observations. 

In the second scenario, the binary could have gone 
through an extra phase of spin-down besides tidal dissi- 
pation. In the sample of Strassmeier et al. (2012) 74% of 
the rapidly-rotating active binary stars are synchronized 
and in circular orbits, but 26% (61 systems) are rotating 
asynchronously of which half have P rot >P r6, mostly giant 
stars. It is suggested that a magnetic wind could have ap- 
plied a braking torque on the stars on the main-sequence, 
the magnitude of which may depend on stellar mass, interior 
structure, and activity. This could indicate that the primary 
and (fully convective) secondary star of WTS 19g-4-02069 
may have been braked with different rates, which could be 
tested by comparing the rotation periods of both stars. Al- 
ternatively, the subsynchronous rotation may have been es- 
tablished during the pre-main sequence phase of the binary, 
through magnetic interaction with a circumbinary disk with 
a central hole (e.g. Casey et al. 1993). Here, the magnetic 
field of the primary could have coupled to the (slower) disk 
motion, slowing down rotation due to angular momentum 
transfer to the disk. Arguably, this would require WTS 19g- 
4-02069 to be young, because any asynchronous rotation 
could be rapidly dissipated through tidal interaction on the 
main- sequence, however we currently have no observational 
evidence to suggest a young system. 



6 CONCLUSION 

In this work we present the discovery of a highly unequal- 
mass eclipsing M-dwarf binary (q=0.27) with P=2.44 d, 
using the high-precision near-infrared time-series of the 
WFCAM Transit Survey. We find stellar masses Mi =0.53 
(±0.02)M Q and M 2 =0.143 (±0.006)M Q , and radii of 
i?i=0.51 (±0.01)i? o and i? 2 =0.174 (±0.006)i? o . The com- 
panion star is therefore in a sparsely sampled and important 



M-dwarf mass-regime for studies of Earth-like exoplanets 
which require accurate calibration of their host star radius 
and mass. We suggest that the low mass-ratio of our binary 
may be explained by the different accretion physics such 
as recently proposed by 2-D warm grid-based (Ochi et al. 
2005, de Val-Borro et al. 2011) or magnetic field braking 
simulations (Zhao & Li 2012), which suggest that short pe- 
riod low-q pairs may be abundant around primaries within 
a large mass range. Alternatively, the binary may have been 
isolated from the birth environment early on through ejec- 
tion or rapid gas removal. Since both stars share the same 
metallicity and age and straddle the theoretical dividing line 
between fully and partially convective atmospheres, a com- 
parison can be made to model stellar atmospheres with the 
same isochrone over a wide span of masses. We find that 
both stars have slightly inflated radii compared to 1 Gyr 
model predictions for their mass, but we argue that future 
work will be required to quantify the effects of star spots 
on the light curve solution. The effective temperature of the 
secondary is consistent with theoretical models, but for the 
primary it is lower by ~350 K, which could be explained by 
magnetic inhibition theory (e.g. Chabrier et al. 2007) if the 
primary is an active star and has a significant magnetic field, 
whereas the secondary may not. The detection of a 2.56 d 
~2% signal in the WFCAM light curve is attributed to sub- 
synchronous rotation of a relatively stable star-spot complex 
at high latitude on the magnetically active primary, suggest- 
ing that its dynamo is more like solar-type stars, with the 
Coriolis force driving up star spots to polar latitudes. 
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s 


ngle 






GJ380 


(35) 


0.670(0.033) 


single 


0.605(0.02) 


single 








s 


ngle 






RRCae-B 


(36) 


WD 


0.1825(0.0139) 


WD 


0.209(0.0143) 




3100(100) 




WD 






NNScr-B 


(37) 


SdB 


0.111(0.004) 


SdB 


0.141(0.002) 








SdB 






GKVir 


(38) 


WD 


0.116(0.003) 


WD 


0.155(0.003) 








WD 






SDSSJ1212 0123 


(39) 


WD 


0.273(0.002) 


WD 


0.306(0.007) 








WD 






SDSSJ1210+3347 


(39) 


WD 


0.158(0.006) 


WD 


0.200(0.003) 








WD 







Table 6. Literature values for the M-dwarf systems used in Figure 151 l9l and [TOl References: (1) Maceroni & Montalban (2004), (2) Becker 
et al. (2008), (3) Dimitrov & Kjurkchieva (2010), (4) Hebb et al. (2006), (5) Cakirh & Ibanoglu (2010), (6) Lopez-Morales et al. (2006), 
(7) Ribas (2003), (8) Irwin et al. (2009), (9) Blake et al. (2008), (10) Vida et al. (2009), (11) Helminiak et al. (2012), (12) Vaccaro et al. 
(2007), (13) Devor et al. (2008), (14) Morales et al. (2009), (15) Irwin et al. (2011), (16) Hartman et al. (2011), (17) Birkby et al. (2012), 
(18) Irwin et al. (2007), (19) Carter et al. (2011), (20) Kraus et al. (2011), (21) Davenport et al. (2012), (22) Creevey et al. (2005), (23) 
Shkolnik et al. (2010), (24) Young et al. (2006), (25) Huelamo et al. (2009), (26) Torres & Ribas (2002), (27) Lopez-Morales & Ribas 
(2005), (28) Beatty (2007), (29) Bouchy et al. (2005), (30) Pont et al. (2006), (31) Pont et al. (2005), (32) Bouchy et al. (2005), (33) 
Fernandez et al. (2009), (34) Tal-Or et al. (2013), (35) Segransan et al. (2002), (36) Maxted et al. (2007), (37) Parsons et al. (2010), (38) 
Parsons et al. (2012), (39) Pyrzas et al. (2012). 
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